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Abstract 

A g-difference analogue of the fourth Painleve equation is proposed. Its symmetry structure and some 
particular solutions are investigated. 

Introduction 

The importance of six Painleve equations and their solutions in mathematics and mathematical physics is 
widely accepted. The discrete analogue of Painleve equations has been first recognized in 0, |8| where they 
appear as Schlesinger transformations of the Painleve equations. After the discovery of discrete analogue of the 
Painleve property which is now called the singularity confinement property[Q, a lot of second order difference 
equations, including g-difference equations, have been identified as the discrete Painleve equations Jl6|, |^|. By 
investigating the Lax pairs, particular solutions including r functions and bilinear equations, and so on||, it 
has been gradually recognized that discrete Painleve equations admit similar properties to original Painleve 
equations. 

The purpose of this paper is to introduce a g-difference analogue of the fourth Painleve equation Pry- 
This equation, called the g-Painlcve IV equation qPw below, shares many characteristic properties with the 
original Pry- In particular, it admits the action of the affine Weyl group of type as a group of Backlund 
transformations. Furthermore, it has an analogue of classical solutions expressible by the continuous g-Hermite- 
Weber functions and rational solutions corresponding those of Pry studied in [jl5| [l0|, |l3) . 

The plan of this paper is as follows. We introduce our qPrv in Section 1 and describe its symmetry structure 
in terms of the affine Weyl group W(A^)- In Section 2, we construct a g-analogue of classical solutions 
along each reflection hyperplane from a seed solution which is described by a discrete Riccati equation. Each 
particular solutions of this class is expressed in terms of a Toeplitz type determinant of continuous q-Hermite- 
Weber functions. In Section 3, we discuss some related topics including the relationship between our gPrv and 
Sakai's Mul.6 system jlTj] and the ultra-discrete limit of gPiv- 



1 Symmetric form of a c/-Painleve equation qP 



1.1 g-Painleve equation gPiv 



In this paper, we consider the following discrete system including three dependent variables / , /i, f 2 and three 
parameters a , a 1; a 2 : 

ao = a , W = ax, 02 = 02, 
/o = a a\fx- 

i ~l~ «0/0 «0«l/0/l 

(1) 



1 - 


F- a 2 f 2 ~ 


F- a 2 a /2/0 


1 - 


H ao/o - 


F- a ai/ /i 


1 - 


F- ao/o - 


F- a ai/ /i 


1 - 


H ai/i - 


F- aia 2 /i/ 2 


1 - 


F- ai/i - 


F- aia 2 /i/ 2 


1 - 


1- a 2 /2 - 


F- a 2 a / 2 / 



/1 = aia 2 / 2 
/ 2 = a 2 a /o 

where — stands for the discrete time evolution. We will also use the notation t(x) = a; when we regard x 
as a transformation of variables. The inverse transformation of t is given by 

= do, ai — ai, a 2 = a 2 , 
/ 2 a ai + a /i + /0/1 



fo = 

k = 
h = 



a ai a 2 a + a 2 / + / 2 / ' 

/o aia2 + ai/2 + /1/2 (2) 
aia 2 a ai + a /i + /0/1 ' 

/1 a2ao + 02/0 + /2/0 
a 2 a aia 2 + a x f 2 + /i/ 2 



We introduce a constant q by setting aoaia 2 = g. Noting that the product /0/1/2 can be regarded as the 
independent variable, we introduce a variable c such that /0/1/2 = ac 2 and c = qc. If we regard fj (j = 0, 1, 2) 
as functions in c, eq. ([!]) thus represents a system of g-difference equations for the unknown functions fj = fj (c) 
(j = 0, 1, 2) with parameters aj (j = 0, 1, 2) such that aoaia 2 = q. As we will see in Subsection 1.3 below, this 
g-difference system has a continuous limit to the symmetric form of the fourth Painleve equation F\y. For this 
reason, we refer to eq.(bj) as the symmetric form of the fourth g-Painleve equation gPiv- 



1.2 Backlund transformations 

The discrete system eq.(|l|) admits the action of the (extended) affine Weyl group W = (s 0: si,s 2 , 7r) of type 

as a group of Backlund transformations. In what follows we denote by W the group generate by the generators 

So, s ii s 2, t and the the fundamental relations 

sj = 1, {s lSt+1 f = 1, 7T 3 = 1, TTSi = S t+1 TT (i = 0, 1, 2), (3) 

where the indices are understood as elements of Z/3Z. 

We define the action of So, s ij s 2> tt on the parameters ao, ai, a 2 by the following formulas: 

so(ao) = a^ 1 , s (ai) = a%ao, s (a 2 ) = a 2 a , 

si(oo) = aoOi, si(ai) = , si(a 2 )=a 2 ai, , , 

S2(ao) = a a 2 , s 2 (ai) = aia 2 , s 2 (a 2 ) = a^ 1 , 

7r(a )=ai, 7r(ai) = a 2 , 7r(a 2 ) = a . 

Note that this definition is a multiplicative version of the standard realization of W on the simple roots aj 
defined through aj — q aj (j — 0, 1, 2). We now define the action of s , sx, s 2 , n on the / variables as follows: 

rt\ t it\ f a o + fo f f v , 1 + ao/o 

so(/o)=Jo, s (/i) = /i— s {f2)=h 

1 + a /o a + /o 

/m f 1 + ai/i . . , . ai + /i 

«i(/o) = /o — Si(/i) = /l, si(/ 2 )=/ 2 



Oi + /l 1 + ai/i (5) 

/m f a 2 + f 2 l + a 2 f 2 

s 2 (/o) = /o^— t-> s 2 (/i) = /i — s 2 (f 2 ) = f 2 , 

1 + a 2 / 2 a 2 + / 2 

tt(/o) = /i, tt(/i) = / 2 , 7T(/ 2 ) - /„. 



We remark that the action of Sj on the a- variables and the /-variables is given by 



Si(aj) = aj a t ZJ , Siifj) = fj ■ 

\ 1 -\~ dij', 



{ij = 0,1,2), 



respectively, where A — (<%■)? -_ ^ s ^ ne generalized Cartan matrix of type A 2 and U = (tiy)-j =0 is an 



(6) 



orientation matrix of the corresponding Dynkin diagram 

A = 



2 -1 -1 
-1 2 -1 
-1 -1 2 



U = 



1 
-1 1 

1 -1 



1 



(7) 



The following theorem can be verified by direct computation. 

Theorem 1.1 The transformations So, Si, S2 and ir of the a-variables and the f -variables, defined by (^), (^), 
generate the extended affine Weyl group W — (sq, Si, s 2 , it) of type A<p . Furthermore, they commute with the 
time evolution t of the fourth q-Painleve equation qPiv- 

By using the action of the extended affine Weyl group W , we can define the Schlesinger transformations T%, 
T 2 , T 3 for gPiv as 

T X = 7TS 2 Sl, T 2 = Sl7TS 2 , T 3 = S 2 SlTT. (8) 

Note that T{Fj = TjTi (i,j — 1, 2, 3) and T\T 2 T 3 = 1. The action of T\ on the variables aj, fj is given explicitly 
as follows: 

Ti(a ) = ga , li(ai) = q~ 1 a x , Ti(a 2 )=a 2 , 



Ti(h) 
Ti{h) 



f 

h 



^ («o + /o)(«o + /o + Qoa 2 / 2 + Q§a 2 /o/ 2 ) 
(1 + aoio)(«o a 2 + a oa 2 /o + /2 + a /o/ 2 ) ' 
1 + ao/o 



(9) 



Ti(f 2 ) = fo 



ao + fo ' 

ala-2 + aoQ 2 /o + / 2 + ao/0/2 
ao + fo + a a a 2 f 2 + a^a 2 / / 2 ' 



The corresponding formulas for T 2 and T3 are obtained by the rotation of indices, since T 2 — irTiir^ 1 and 
T3 = ttT^" 1 . Each of these Schlesinger transformations commutes with the time evolution of gPrv and can be 
regarded as a version of the third q-Painleve equation. For example, eq.(^|) is rewritten equivalently as 



1 + a fo qc 2 1 + a .fo 



ao + fo fofi a + fo ' 

rr-i(f^f ai + h _ qc 2 ai + /i 

J i \Jo — J 2 — j- — -j-r t~. 7~- 

1 + ai/i /0/1 1 + ai/i 



(10) 



This is an analogue of the fact that some discrete Painleve equations arise from Backlund transformations of 
the Painleve equations |, |, §, [§, [l2| . 

1.3 Limit transition to P IV 

2 Icy 

By introducing a small parameter e such that q = e~ e ' , we set 



a« = e- £Q */ 2 , /i = -e-^ (j = 0,1,2). 



Then we have 



ft - fi = e 2 (tpi(ip i+ i - ifi+ 2 ) + a,) + 0(e 3 ) 
for i = 0, 1, 2. Passing to the variables aj and ipj, we define the derivation ' by 

x' = lim — (x — x) 

E^O £ 



(11) 
(12) 

(13) 



for a function x in aj and ipj. Then we have /, — /,; = s 1 Lp' i + 0(e 3 ). Hence we obtain 



a' 


= 0, ai = 


= 0, 


a 2 


<Pa 


= Po(<Pi ~ 


^2) 


+ ao 


<f'i 


= ^1(^2 - 




+ Oil 


•r-2 


= ^2(^3 - 


Vi) 


+ a 2 



(14) 



This differential system is the symmetric form the fourth Painleve equation Piv [[l9| [j], In fact, under the 
normalization tpo + ip\ + ip 2 = t, «o + a i + a 2 = 1, the system ( fL4| ) is equivalent to the second order differential 
equation 

y" - ^(yf + \f - ^ + (£ + ai - a 2 )y - g (15) 

for y = ipa, where ' = d/dt. Through the limiting procedure, the Backhand transformations Sq, s±, s 2 , 7r for our 
qPrv also pass to those for the symmetric form of Piv ( |l4| ) such that 

Si(aj) =a j -a i a ij , Si(<pj) = <p, + , . 

¥>i (z, j =0,1,2) (16) 

as in El P. 



2 Particular solutions for 



IV 



In this section, we investigate g-analogue of classical solutions of Piv and determinant formulas for them. Recall 
that Piv has two classes of classical solutions [|l5| |l0|, [l4| , those of hypergeometric type, expressed in terms of 
Hcrmite- Weber functions and the rational solutions. 

2.1 Continuous g-Hermite- Weber functions as seed solutions 

In order to show the parameter dependence explicitly, we rewrite gPrv (Q) by operating T"T 2 on ([!]) (y, JVsZ). 
We denote T"T 2 N (fi) = fi(c;v,N), and we abbreviate the unnecessary arguments depending on the context. 
Notice that 

T?T 2 N (a ) = a q\ T»T» {a x ) = a iq ^+ N , T»T 2 N \a 2 ) = a 2 q- N , T?T 2 N (c)=c, (17) 

and 

t{ai)= ai , (i = 0, 1,2), t{c)=qc, (18) 

where 

q = a aia 2 , c= ( — ~ — ) • ( 19 ) 
The g-Painleve IV equation qPiy with respect to the variable c is expressed as 

tt \ N f , ^ + a2g~ jV /2(c) +a 2 a og ' / - Ar / 2 (c)/o(c) 

f (qc)=a a 1 q /i(c — — n~FT\TT\ > 20 

1 + ao9 y /o(c) + aoaig iv /o(c)/i(c) 

, / n -fjr^ 1 + a o ol/ /o(c) +a ai q N fo(c)fi{c) 

f 1 (qc)=a 1 a 2 q /2(c)— _^ +jv FTTT f H f c V 21 

ff x , M 1 + aig-" +A 7i(<0 + aia2g- / /i(c)/ 2 (c) 

f 2 {qc)=a 2 a Q q /o(c) — ZWTT\ 1 ^=jv7 rTTTT " 22 

l + a 2 o JV / 2 (c) +a 2 aog^ JV /2(c) /o(c) 

where the arguments ^, AT for /$ are suppressed so that /i(c) = fi(c;v, N). In what follows, we use similar 
abbreviations 

/ t (<Z fc c) = /,(<z fc c;^,7V), A(i/ + A!)=/ 1 (c;i/ + A! > JV), /,(JV + fc) = /„(<:; 1/, JV + fc) (23) 



for k E Z. 

First let us consider the case N = 0. It is possible to specialize the variables as 



/ 2 = -l, o 2 = l, (24) 



consistently. In fact, eq.([22|) becomes trivial by this specialization, and eqs.(|20|) and (|21j) are reduced to a 
discrete Riccati type equation 

h M = -JL- ii-Wff)-««r* (25) 



and /o(c) = ■ By putting / x (c) = |^ in eq.(|25j), we can solve eq.(|2 



as 



g G„(gc 

A c ) = U n u\ ' 26 

dog Cri/(CJ 

with Gj/(c) satisfying the linear g-difference equation 

G,(g 2 c) = (1 - g 2 c 2 )G,(gc) + alq 2v c 2 G v (c). (27) 

Let us derive contiguity relations to be satisfied by G v (c) in the direction of v. For this purpose, we consider 
the Schlesinger transformation T\ (|9|), 

/ (^ + i) = h(v) 



1 + a q»f {v) 

a 2 a 2 q 2v - N + a a 2 q»- N f (is) + f 2 (u) + a q» f {v) f 2 (v) ' 



(28) 



tl l1 n 1 + a q v fo(v) . . 

1) = / 2 (f) - ry-r, 29 

/ a (f + i) = /oM 

X aog , ' + /oM+aoa2g^ A '/2(^+aga 2 g 2 ^/oH^H■ 1 j 



Applying the specialization ( |24| ) and putting TV = 0, cqs.(p8|)-(|30D are reduced to a discrete Riccati type equation 



/ 1 M-q g'-+ 1 c 2 
a^q v j\(v) - qc 2 



h(-+i) = - jn :i; f T _ 2 . (3i) 



Substituting eq.(|26|) into eq.(|3l|), we obtain contiguity relations for G v (c), 

G v+1 (qc) = G v (qc) + a 2 q 2v c 2 G v (c) , (32) 
G u+X (c) = G v {qc) + c 2 G„(c). (33) 

In particular, we obtain a three-term relation in the direction of v, 

G v+2 (c) = (c 2 + 1) G v+1 (c) - c 2 (l - a 2 q 2v )G v (c). (34) 



We note that eqs. (|27j) and ( |34| ) are derived from eqs.(|32|) and 
From the above discussion, we have: 



Proposition 2.1 qPiv admits 1 -parameter family of particular solutions given by 



2 v G ^( c ) t 1 G »(l c ) 



fa = c a q v A = -. , , h = -1 (35) 

G„(gc) a Q q u G v (c) 

for N = 0, where G v (c) is a function satisfying the contiguity relations and 



It is interesting to note that eq. (J27|) admits polynomial solutions in c if ao = q, 

G (c) = 1, 

Gi(c) = c 2 + l, 

G 2 (c) = c 4 + (1 + q 2 )c 2 + 1, 

G 3 (c) = c 6 + (l + < Z 2 +g 4 )c 4 + (l + 9 2 +<zV + l, 

It is not difficult to check that 

-E f2 if''£\ ^ = E □ , e2iy " 2fc (" = 0,1,2,. •■) (36) 
^ {q 2 ;q 2 )k(q 2 ;q 2 )v-k k~o q 

satisfy eqs. (^||) and (|33|), where [£] is the g-binomial coefficient. The generating function for these polynomials 
is given by 

1 _^ G„(c) (XV 

(cX;q 2 )(X/c;q 2 ) ^(q 2 ;q 2 ), W ' 1 ' 



u=0 

The polynomials H v (x) = cr u G v (c), X — £1X6 called the continuous q-Hermite polynomials Jll|, |). In this 

sense, c~"G v {c) for a ^ q l (I 6 Z) may be regarded as a g-difference analogue of the Hermite- Weber functions, 
to which hereafter we will refer as "continuous q- Hermite- Weber functions" . 

We note that it is possible to obtain "higher order" solutions for N 6 Z by successive applications of T2 
on the solutions obtained aboveQ- These solutions are expressed rationally in the continuous g-Hermite- Weber 
functions. 

From the rotational symmetry 7r of qPw, it is possible to apply other specialization, (an,/o) = (X, — 1), 



v = or (<2i, /1) = (1, —1), v — N = on gPiy(|20|)-(|22[) and perform the same procedure as discussed above. 
Therefore, we obtain 

Theorem 2.2 When <n — q k for some i = 0,1,2 and k € Z, the fourth q-Painleve equation qPiv ^) admits 
a 1 -parameter family of particular solutions which are expressed rationally by the continuous q-Hermite- Weber 
functions. 

Explicit description of these solutions will be given in the next section. 

Note that the hypersurfaces a,; = q k (i = 0, 1, 2; k 6 Z) in the parameter space correspond to the reflection 
hyperplanes of the affine Weyl group W = W{j^). This is an analogue of a well-known result by Okamoto 
for the classical solutions of the continuous Piv 

2.2 Toeplitz determinants 

2.2.1 Determinant formula and bilinear equations 

In this section, we prove the following determinant formula for the particular solutions obtained in the previous 
section. 



Theorem 2.3 Let G v {c) is a solution of eqs. and \3!y. For each N e Z>o, we define an N x N Toeplitz 
determinant 

r N (c)=det(G„- l+3 (c)) ij=1 ... N (NeZ >0 ), $(c) = l. (38) 

Then, 



v 2 0at + i(c)0at(?c) 



<l>N +1 (qc)4>N(c) 

h{c;v,N) = -— — — — — — - , (40) 

a q" + " (f> l/ N+1 {c)(f> N (qc) 



^ N {qc)^ N ~\c) 



f 2 (c;u,N) = -q N ™ , ^lu-iV., («) 



lr The situation is a little delicate for N < 0, but we will discuss in Section 2.2.3. 



satisfy qP Iv ^2^-^i) with a 2 = 1. 



We note that the case of N = agrees with Proposition 2.1, 

Theorem 2.3 is a direct consequence of the following proposition: 

Proposition 2.4 <j) N (c) satisfies the following bilinear difference equations, 

alc 2 q 2 »<P N {<lc)<P N+1 {c) + 4>n(cWn+i(Qc) = q-^^+M^N 1 ^), 
aW^qcW+^ic) - ^(cW+^qc) = (a 2 q 2 » q 2N )d> N+l {qc)<t> v +\c), 
#&(?<0^(c) - q 2N ^ N + M^N(lc) = c 2 (alq 2 » - q 2N ) ^ N +i (c)^ +1 (qc), 
r N +1 (qcW N (c) ~ <ftf X {c)<f> N {qc) = -^(gc^+^c), 

r N +1 (qcW N (c) ~ q 2N r N +1 (cW N (qc) = -T^^+M^N-M 

a Q q c 

q 2N r N (cW N+1 (c) - ^(0)^(0) = -a 2 c 2 q 2 ^rN + i(Q- X cW N (qc), 

Pn(cWU c ) + c2 ^ 2u ~ q 2N W N+ i(cW N +1 (c) = q-^Mt-^WMqc), 

a 2 Q q 2 »^ N (cW N +1 (c) + r N + +^ N -i(c) = alq^^+'iq-'cW^qc). 



(42) 
(43) 
(44) 
(45) 

(46) 

(47) 
(48) 
(49) 



Proof of Theorem We show that fc(c) {i = 0,1,2) defined by eqs.(]3J)-(|4lJ) satisfy eq. (|20|) . Substituting 
eqs.(p9)-(41) into the numerator of right hand side of eq.(pfl), we have 



• f 2 (c) + a q / 2 (c)/ (c) = 1 + g / 2 (c) l + a g c 



^ +1 (gc)^(c) 



/2(C) 



0W?c)#v(c) + «^^c 2 ^ +1 (c)^((zc) 



^+i(?c)^(c) 



= l + q- N f 2 (c)q 



<t> N+ i(qc)cf> u N (c) ' ' (f N {qc)4> N+1 {qc) 



4> N (l0)r N+ l(qc) a ° Cq r N {qc)r N+1 {qcY 

where we used bilinear equations (E3) and fli?]). The denominator is calculated by using eq.(|43|) and 



as 



l + a q v f {c)+q N+1 f {c)h{c) = 



q- 2N ^ N \c) 



'N+1 



(q 2 c) 



4> N+1 {qc)<t> N 1 {q°) 



(50) 



Then, the right hand side of eq. ([20|) reduces to 

l + q- N f2(c)+a q»- N f 2 (c)f (c) 



JV+1 



h{c)- 



l + am v h{c)+q N+1 h{c)fi{c) 

v-l 



q 



N+1 



1 



a q 



u+N-l iv 



^N+i(lcW N (c) 2 2 a^-jy) </>^ +1 (c)^(g 2 c) 2N 4> N + 1 (qc)4> N (qc) 

r N +i(cWN\qc) ' a ° c q r N {qcw N+1 (qc) ' q 



^ _1 (c)^ +1 (9 2 c) 



2 „+2 ^Ar+l(gc)^Af(g C) 

fl o c <? -n; — TTsmi/ — = /o(<7 c ). 



r N +i(q 2 cW N (qc) 



is checked by using the bilinear equations (|43|), (47), ( |45| ) 



Thus we have shown that eq. (|20|) follows. Eq. 
and(^). Eq.(p2|) follows automatically. 

Remark. It should be noted that both 1 + ctifi and 1 + /,-/oj (i = 0, 1,2) admit multiplicative formula with 
respect to 4>. In fact, we have, 



1 + ao9 I/ /o(c; v,N) = q 



-2N 



<l>N +1 (qcW N (c) ' 



(51) 



1 + aiq 



-1/+JV 



/i(c; v,N) 



,2 g 2(,-l) _ g 2iV ^-l l(gc )^( c ) 

<^ +1 (c)0^V)' 



l + a 2 ^ J V 2 (c; i/, AO 



ag g 2(,-l) 



' ^ _1 (c)^(gc) 



and 



r N+1 (gc)r N (c) ' 



ao<T 
/ 2 (c; t/.JV) 

-JV 



2, >.\ 2J2(v-l)s ^jV+lXgMjvjgg) 

^+i(c)0jv (?c) 



1 



$v+i(g c )<^-i( c ) 



which can be verified directly by using eqs.(|39|)-(f4l|) and the bilinear equations 



(52) 
(53) 

(54) 
(55) 
(56) 



2.2.2 Proof of Proposition 2.4 



Our basic idea for proving Proposition 2.4 is as follows. Bilinear difference equations are derived from the 
Pliicker relations, which are quadratic identities among determinants whose columns are shifted. Therefore, we 
first construct such "difference formulas" that relate "shifted determinants" and <^(c), by using the contiguity 
relations of G„(c). We then derive bilinear difference equations with the aid of difference formulas from proper 
Pliicker relations. We take eq.([42|) as an example to show this procedure explicitly. For other equations, see 
Appendix A. 



Let us introduce notations, 



N 



G„(q m c) = G™ 

Gm 

Grn 
' ' ' 



v+N-l 
m 

u-N-2 



Gn 
„ 



J v-N+1 

where k m denotes a column vector, 



G 



u-N+2 



Gn 



i fn x m 



N- V 



(57) 
(58) 



/ G 

G 



u+k 
i 

v+k-1 



(59) 



\ (-im 

\ ^v+k-N+l 

Here the height of the column is N , but we use the same symbol for determinants with different size, since there 
is no possibility of confusion. 

We next construct a difference formula. 



Lemma 2.5 (Difference Formula I) 



| m+i/ j o ro , l m ,---, N-3 m , N-2 
| 1 m+l/ ) Q m ^..^ N _ 3 ™ N _ 2 
where k m/ is a column vector, 

( 



-2(iV-l) 2(N-1) -2(JV-l)(m+y-l) ±y,m+l 



J N 



mi _ -2(iV-l) 2(iV-l) _-2(AT-l)(m+i/-l) jv,m+l 



1 U u+k-l 



(60) 
(61) 

(62) 



+k-N+l 



Proof of Lemma 2.5. We use the contiguity relation j3^), which is rewritten as 

G™ +1 + a 2 ,<? 2(l/+m) c 2 G™. 



Gm+l 



(63) 



in the present notation. Then we have 



ii/,m+l 



which is eq. 
get 



Gm+l 
v-X 



✓"fm+1 
£fm+l 

Gm+l 
i/-l 



Gm+l 
I/+1 



m+l 



G 



W+l 

/""Vm+l 



^iz+AT-l 

✓^m+1 

U y-iV-2 



y^mi+1 



G™ +1 
G 



rn+l 



^zy+iV-l ~ U zy+iV-2 

y^-f ??i+l ✓~*m+l 

°V+iV-2 ~ U zy+iV-3 



m+l 



W-JV+1 



G 



m+l 



Gm+l 
v-1 



m+l 
m+l 



a 2 c2(? 2(^+m) G m 
a 2 c 2 q 2(u-l+m) G m_ x 



n 2 n 2{y+N-2+va)„2 r <m 
a a c W+iV-2 

n 2 n 2{y+N-S,+va) r 2 r <m 



W-jv+i 



„2„2„2(iy-AT+l+ m ) r( m+l _2,.2(j/-l+m)„2/'-fm 

a c q <-V-jv+i a o1 C W-l 

/^m+1 /~fm /^f 

^2^-tm+l .r 



2(JV-1) 2(AT-l)„2(iV-l)( I /+m-l) 



T i/-1 



G 



v+N-2 
m 

v+N-3 



n 2{N-l) r m 



a 2(JV-l) c 2(W-l)„2(AT-l)(i/+m-l) | n m+l/ n m 



'q--- — ' |0 m+1 ', m , l m ,---, N-3 m , N-2 m |, 
In the second line of the above calculation, adding the second column to the first column, we 



< i/,m+l 
9 JV 



/^m+1 



/"fm+1 /^(m+1 sim+l 



^im+1 



U ^+Af-2 



W-JV-2 — W+N-3 



/-im+l 
<J i/-iV+2 



/"fm+1 
W-Ar+2 



^V-iV+l 



G 



m+l 



Gm+l 
z^-1 



from which we obtain eq. (|6l|) by using eq.(|63] 
Wc then consider the Pliickcr relation, 



(64) 



= l m+1/ , l m ,---,N-l m | x |0 m , l m , 

- \ip u m , l m ,---,N-l m | x |l m+1/ , 1™, 
+ |l m+1 ', m , l m ,---,N-l m | x l m , 



•,N-l m , N 1 
,N-l m , 1ST 
•,N-1 TO , N 1 



(65) 



where 



/ 1 \ 





V o / 



By expansion with respect to the column <pi, the Pliicker relation (65) is rewritten as 



(66) 



= 



q 2 \0 m+v , m ,-- 
|-l m , m ,---,N 

I -i 777+1/ nm -i m 
M- J } J- 3 * 



,N-2 m , N-l m x |0 m , l m ,---,N-l 



N 



x 1 



m+l/ 



,N- 



,N- l m x |0 r ' 



-,N- 2 m , N-l' 



where we have used 



m+lV -i rn 



<pi, 1 

q 2 G™ +1 



i™ ... n - r 



9 



2gm+l 



G 



v+N-1 
m 

v+N-2 



-iV+1 °V-JV+2 



n n 2N/~<m+ 



GI: 



G 



„2N rrm+l 

1 °V-JV+1 °V-A»"+2 



iz+AT-2 



G!; 



^2 Q™ 



N 



Then, we obtain by using Lemma 2.5 



= fl 2 X flQ 2(A r -l) c -2(Ar-l) (Z -2( J V-l)(, +m -l) x 

jLiz-l.m -2iV -2iV -2iV(i/+m) il/+l,m+l 
X Uq C Q ' ' 



iV 



'0 



-2JV„-2iV_-2iV(i/+m-l) 



iV ' 



(67) 



which is the same as eq.(|42j). This completes the proof of the bilinear equation (|4 
2.2.3 Determinant formula for negative N 



Proposition 2.3 is a determinant formula for the solutions with N € Z>o, which are obtained by successive 
application of T 2 on the seed solution described in Proposition 2.1. In order to obtain solutions for N € Z<o, 
we have to apply T 2 _1 on the seed solution, however, we find that this procedure collapses. Therefore, we 
construct "seed and higher solutions" for negative N so that bilinear equations described in Proposition 2.4 
hold for all N 6 Z<o- To find the seed solutions for N = — 1, we put N = — 1, <^i(c) = G„(c) and <Pq(c) = 1 
in eqs.(El)-(E^), we obtain a set of contiguity relations to be satisfied by G„(c), 



a 2 c 2 q 2v G v (qc) + G v {c) = q 2 G „_i(c), 
G v {qc) = G v+1 {c) + c*G v+ x{qc). 



(68) 
(69) 



Similarly to the case of N > 0, the following three-term relations in the direction of c or v are derived from 
eqs.© and ©, 



G„(c) 



1 



1 



c 2 )GAqc) + a 2 q 2 »c 2 G„(q 2 c), 

2 2v\ 



G v {c) = - (1 + c 2 ) G„ +1 (c) - - (1 - a 2 ,? 2 ") G„ +2 (c) 



(70) 
(71) 



Theorem 2.6 Lei G v (c) is a function satisfying eqs. ^6dj ) and IpQ). For each —M = N G Z<o, we define an 
M X M Toeplitz determinant 



Then we obtain a determinant formula for N < as follows. 



<^(c) = det (G„_ i+j (c)) ... (-M = N e Z <0 ), 



(72) 



Then, <j) v N {c) satisfies the bilinear equations 

02 = 1. 



and therefore eqs.^^-j^ safe/?/ qP Iv ^2dj)-^) with 



Since this theorem is proved by similar procedure to Theorem 2.3, we omit the details. 
Remark If we parametrize as ao = o M , 

G„(c; q) = c 2 q- 2u G_ {l ,_ 2ll) (c; 1/q), 

satisfies eqs.(p§) and (p9) formally. 



(73) 



2.3 g-Okamoto polynomials 

In the previous sections we have discussed the particular solutions on the reflection hyperplane of the affinc 
Weyl group W = W(A^) in the parameter space. It is easy to see that qPiy (|lj) with normalization condition 
( |l9| ) admits a particular solution given by 

(fo, fi, f 2; ao, 0,1,0.2) = (xp^ 1 ,xp^ 1 ,xp^ 1 ;p^ 1 ,p^ 1 ,p^ 1 ), x = c 2/3 , p^q- 1 ? 3 . (74) 

Note that this parameter corresponds to the fixed point of 7r, namely Dynkin diagram automorphism of the 
affine Weyl group W(A^). By applying the Backlund transformations on this seed solution, we get a series of 
rational solutions on the barycenters of the Weyl chamber of the affine Weyl group W(A^). 

Theorem 2.7 For the parameters 

T™T?(a , ai ,a 2 ) = (p- 3 ^ 1 , p 3 ^ 3 ^ 1 , p 3 ^ 1 ) , m,neZ, (75) 
we have the following rational solutions of qPiv equation 

l — 2n—l Qm-\-l,n 

(xp 2 )Q m +l,n+l{x) 



T{ n T?(f ) = X P " 



Qm+l,n{x)Q m +l,n+l{xp 2 ) ' 



rpmrpml f \ _ --.m+n-1 Q"i+l,n+l(xp )Q m ,n{x) , . 

1 i 1 2 \Ji)—xp — TYFi 1 ^2T' v . lo ) 

Qm+l.n+l{x)Q m .n(Xp Z ) 
rpmrpn^j^ _ X pn-2m-l Q m , n ( X P 2 )Qm+l,n(x) 

(x)Q m+1 , n (xp 2 ) 

Here Q m ,n{x) are monic polynomials defined through the recurrence relations, 

Q m -l,n(x)Qm+l,n+l(xP~ 2 ) = <3m,n(l)<3m,n+l(^ 2 ) + Xp~ 2m ~~ 2n Q m , n+ i{x)Q m , n {xp~ 2 ), 

Q m +lAi(x)Qm,n+l(xp~ 2 ) = Q m +l,n+l(x)Q m ,n( X P~ 2 ) + Xp 4m ~ 2n ~ 2 Q m ,n (x)Q m +l,n+l (xp~ 2 ) , 
Qrn+l,n+l{x)Q m ,n-l(xp~ 2 ) = Q m ,n (x)Qm+l.n (xp~ 2 ) + Xp 4n ~ 2m ~ 2 Q m+1 ^ n (x)Q m . n (xp~ 2 ) , 

with initial conditions Qo.o(x) = Qi^(x) = Q\ t i{x) = 1. 

We call the polynomials Q m , n {x) the q-Okamoto polynomials. Some examples are as follows. 
<22,oO) = x 2 + {p 2 + l)x + 1, 

Qs,o( x ) = x 6 + x 5 (p 6 +p 4 +p 2 + 2+p- 2 ) + x' l (2p 8 + 2p 6 + 3p 4 + 3p 2 + 3+p- 2 +p- 4 ) 
+x 3 (p w + 2p s + 3p 6 + 4p 4 + 4p 2 + 3 + 2p- 2 + p- 4 ) 

+x 2 {2p 8 + 2p e + 3p 4 + 3p 2 + 3 + p- 2 + p- 4 ) + x(p 6 +p 4 +p 2 + 2+p- 2 ) + l, 

Q2,l(x) = X + l, 

Q 3 ,i(x) = x 4 + x 3 (p 4 +p 2 + l+p- 2 ) + x 2 (2p 4 +p 2 + 2+p- 2 )+x{p 4 +p 2 + l+p- 2 ) + l. 

Similarly to the conti nuou s case, g-Okamoto polynomials admit determinant formula of Jacobi-Trudi type. The 
proof of the Theorem 2J as well as the Jacobi-Trudi type formulas will be given in our next paper 0, where we 
will discuss the r functions in a more general setting. 



3 Discussion 



3.1 W(A^) x VT(A^) symmetry and comparison with Sakai's Mul.6 system 

Along with the transformations Si and it of W^A^ 1 ), we define the transformations wq, w\ and r acting on 
C(a i; fi (i = 0, 1, 2)) as follows: 

a a 1 (a 2 a + a 2 f + / 2 /o) , , , 1 + o /o + a-oai/o/i 

w o(/o) = —T7 : 1 ; , 7 x . ^i(/o) = 



/ 2 (a ai + a /i + /o/i) ' a ai/i(l + a 2 / 2 + a 2 a / 2 / ) 

aia 2 (a ai + a /i + /o/i) , , , 1 + &i/i + »ia 2 /i/ 2 

«Wi) = — 77 ; 1 rTTT ' w i(/i) ~ 



/o(ai0 2 + ai/ 2 + /i/ 2 ) ' aia 2 / 2 (l + a / + a ai/ /i) 

, , \ a 2 a (aia 2 + ai/ 2 + /i/ 2 ) , , . 1 + a 2 f 2 + a 2 a / 2 /o 

^o(j2j = —77 ; r -77r < ">i(/ 2 ) = 



(77) 



/l(o 2 a + a 2 / + / 2 /o) ' 0200/0(1 + 01/1 + 0102/1/2) 

= 7-, wo(oj) = toi (oi) = r(ai) = Oj, (i = 0,1,2). 

Lemma 3.1 T/ie transformations wq, Wx and r generate the extended affine Weyl group W(A± ), namely we 
have 

Wq = w\ = r 2 = 1, rwo = w\r. (78) 



Moreover, this action ofW(A^) = (wo,w±,r) commutes with that ofW(A%) — (so,Si,s 2 , 

Note that the discrete time evolution of the qPiv system is a translation of the W(A^ ), that is, 

t = T i = rw . (79) 

Lemma 3.2 The representation ofW(A^) x W(.Ai ) is equivalent to Sakai's system Mul.6. 

Proof. Let x,y,z and do, &i, a 2 , 61, 60 = dodxd 2 /bx be Sakai's homogeneous variables and the parameters |jl7|. 
Then his representation of W(A^) x W(A^) is given as follows: 

7T = (T 4 , S 2 = irSiTT 2 , Sq — 7T 2 S 1 7T, V = <7 3 , = rW\ T. (80) 

Where <r, s\ and 101 is given by, ^ 

cr(a ) = ai, a (ax) = a 2 , er(a 2 ) = a , <r(b ) = &i, cr(&i) = 6 , 
cr(x) = d 2 xy(z — a;), cr(y) = ~bxyz(x + y — z), a(z) — a 2 x(x — z) 2 , 



sx(a ) = a ax, s 1 (d 1 ) = d 1 1 , sx(a 2 ) = d 2 dx, si(ftj) = (i = 0, 1) 
si(x)=x, sx(y) = axy, sx(z) = &xz, 

w 1 (d l ) = di, (i = 0,1,2) wi(6o) = bob?, iUi(&i) = frf 1 , 

= + a 2 x), wi(y) = -biy(y + a 2 a;), k;i(z) = 2r(a 2 ar - &iy). 



(81) 



Introduce variables ft (i — 0, 1, 2) as 



f - x f - a x ~ z f - biyz (*?\ 

71 — , J2 - «2 , Jo—- — -, r. (oZ) 

z y a 2 x(x — z) 



Then we have 



o-(/o) = T , o-(fx) = -=-, <r(f 2 ) = T , 
7i 72 Jo 

,?v ;Oi(l + /i) fx ,7. Oi + / 2 

si(/o) = /o— : j — , si(/i) = — , si(/ 2 ) = (83) 



ai + fx oi ' " 1 + f- 



ff s 1 + /o + /o/i ,7. 1 + fi + /1/2 I + /2 + /2/0 

w U/o) = -~t: — ? ttt w h/i) = tt: — ; r^7 ^1(72) ~ 



/i(l + / 2 + / 2 /o) / 2 (l + /o + /o/i) ' /o(l + /i+/i/ 2 ) 

This representation is equivalent to our representation by the relation /j = Oj/j, 5.^ = a?. I 

2 The generators cr, si, uii here correspond with o"( 12 345o) 1 w i an< i w 'x 111 jl^j • ^ n error i n the formula of cr( 12 345o) is corrected. 



3.2 Ultra-discretization of gPiv 

In Section 1, we have mentioned that <?Piv admit a continuous limit to the symmetric form of Prv- There is 
another interesting limit, which is known as "ultra-discrete limit" [fl8|. We put 



a t = e Al/e 



i = 0,1,2 



and take the limit e 



-0. By using the formula, 
lim elog ( e^ + 



max(A, B, 



(84) 



gPrv(M yields 



F =t(F ) = Aq + A-l +Fi + max(0, A 2 + F 2 , A + A 2 + F + F 2 ) 

- max (0, A Q + F , A x + A + Fi + F ) , 
F 1 =t{F 1 ) = Ai + A 2 +F 2 + max(0, A + F , A ± + A + F ± + F ) 

- max(0, Ai+Fi, A 2 + A x + F 2 + Fi), 
F 2 = t(F 2 ) = A 2 + A Q + F + max (0, A 1 + F lt A 2 + A Y + F 2 + Ft) 

- max (0, A 2 + F 2 , A + A 2 + F + F 2 ) , 

Ai = t(Ai) = Ai, 



(85) 



which we call the fourth ultra-discrete Painleve equation(uPrv)- Simultaneously, Backlund transformation 
is ultra-discretized as 



s^Fj) = Fj + Uij (max(Ai, Fi) - max(0, A. t + F,)) , 

s i (A j ) = A j -a ij A i , ir(X l ) = X l+1 (i = 0,1,2), X = F,A 

where A = (aij)i,j=o,i,2 and U — (u,j)i J= o.i,2 are given by eq.(Q). Then we can verify the following: 



(86) 



Proposition 3.3 The transformations sq, s\, s 2 and n of the A-variables and the F -variables, defined by 
eg.([$^) generate the extended affine Weyl group W = (sq, s%, s 2 , 7r) of type A^ ■ Furthermore, they commute 
with the time evolution t of the fourth ultra- discrete Painleve equation uPjy. 

A Difference Formulas and Pliicker Relations 



In this appendix, we provide with data which arc necessary for the proof of Proposition 2.4 



We first note that it is possible to express 4>^ in Section 2.2 as Casorati determinant with respect to m as 
follows. 



Lemma A.l 



N 



is rewritten as 

Gm 
v 

Gm+1 
„ 



Gm 

Gm+1 



G 



m+N-l 



G 



m+N—X 
v-X 



N-X 



iv.ra 



n 



2_2(n-fc) 



\ \alq 2( - n -V - 1 



N-k 



rim 

^v-N+X 
rim-\-X 



rtm+N — 1 
^*v-N+X 

Gm 
v 

Gm— 1 
„ 



G 



m-N+X 



rim+1 

Gm 
„ 



G 



m-N+2 



rim+N-X 
rnn+N-2 



(87) 



Proof of Lemma A.l. We use the contiguity relation (33), which is rewritten with the notation introduced in 
eq.(|57]) as 

G™ +1 = G™ +1 - q 2m c 2 G™. (89) 



For k = 2 to N, subtracting (j — l)-st column multiplied by q 2m c 2 from j-th column of eq.(|5g|) for j — N, ■ ■ ■ k, 
we obtain 



G 
G 



G 



G 



v-N+1 

G™ 

Gm 
v-\ 

n 

v-N+l 



fim n 2 

Gm „2m 
v ~ 1 



G 



u-N+2 

Gm+1 
v 

Gm+1 
., — 1 



^iz-jV+l 



T v-\ 



L ^iz-jV+l 

v"~Tm+l 

f~im+l 



G 



m+1 



G 

Gn 
„ 



v+N-1 

-N-2 



1 C 

q c i^r v+N _ 3 



J1 2m„2r'm 



which yields eq.(E 



G'; 



Gm+1 
v 

Gm+1 



/~im-\-l 



r<m+N-l 



^tm+N-1 



(90) 



T u-N+1 

by taking the transposition. In order to derive eq. (|8S|) , we use the contiguity relation. 



Gri 



1^2(^-1) _ I 



1 



Gn 



(91) 



which is obtained from eqs.(|32|) and (|33|). For k = 2 to N, subtracting (i — l)-st row multiplied by a 2 g 2(„ 



from j-th row ofeq.@ tor i = N,---k, we obtain 

G™ 



! /A/J7 

»7V 



G 



G 



J/-AT+1 



I/-JV+2 

a 2 2(„-iV + l)_ 1 



"^iz-JV+l 



a 2 g 2(,-iV+l) 

a 2 g2 (*-iv+i) _ 



JV-l 

n 



a 2q,2(*/-;) _ 



Gi: 



iz-JV+2 
m— 1 
v-N+2 



r m+X 
a 2 2(„-lv + l)_ 1 



m+1 



G 



/^fm+l 



G 



Gm 
v 

Gm— 1 



rim 

^v-N+2 
C V-jV+2 



v-N+2 

Gm+1 
v 

Gm 



f~im-\-\ 
U v-N+2 



G 



v-N+2 



m+N-1 



G 



/-tm+N-l 
L, v-N+2 
f~im+N-2 
U v-N+2 



Qm + N-1 
Qm+N-2 



^m+W-l 
U i/-iV+2 
^im+N-2 
U u-N+2 



Qm+N-1 



(~im+N-l 
U i>-Ar+2 



s^m+N-l 
W-JV+1 



i + N-\ 

a'iq 2 ^- N + 1 )-l 



N-l 

n 



a 2 q 2 ^ 



Af-i 



Gm— 1 



Gm+1 
v 

Gm 



rtm-N+2 r<m-N+3 

Gm—N+l r^m — N 



rirn+N-2 

Gm — 1 
v 

Gm 



which is eq.(|S8|). 



I 



In view of this lemma, we introduce the following notations. 



iv.rci 

9n = 


G™ 

rim 


W+l 

G™ •• 


/~im 
' ^v-N-2 


<t>N = 


rim 

G™ 


( - T i/-Af+2 ' ' 


■ G™ 


i i/.m 

= 




s-im+N-l 
/~im+l 

G™ 


s~im+N-l 
rrm+N-l 




fim-N+l 


fim-N+2 


G? 



= |0 m , l m , 



-»m— 1 



• • , n - r 



■ m— 1 



-N 



= |0„, 1 V ,---,N-2 V , N-l„| 



where k m , k m and are column vectors given by 



k m = 



( G 

G 



v+k 
i 

v+k-1 



\ 



k™ = 



V G™ +k _ N+1 / 
respectively. We also use auxiliary column vectors 



r~vm-\-N — 1 
\ U v+k 



\ 
I 



k„ = 



/ G™ +fc \ 



( 



G 



v+k 



q 2 G™ 



\ 



( 



+fe-i 



k mt 



G 
9 <J ,. 



f+fc 

'm-f 

f+fc 



/ G™+ fc 

Jlfim+k — 1 
^ g2(JV-l)grm+/c-AT+l 



\ 1 ^v+k-N+l ) 

Now we give difference formulas. 
Difference Formula I 

| m+l/ } Qm N-3" 1 , N-2" l | =a - 2 ( JV - 1 ) c -2(JV-l) 9 -2(JV-l)(m+ I ,-l) 

\l m + v , m , l m ,---, N-3 m , N-2 m | =a - 2 ( Ar - 1 ) c -2(w-i) ? -2(w-i)(m+ I ,-i) 
Difference Formula II 

|0 m , m+1 , l m+1 ,---, N-3 m+1 , N-2 m+1 | = <^"\ 
|l m , m+1 , l m+1 ,---, N-3 m+1 , N-2 m+1 | = cV m <^ m . 
Difference formula III 

,---, -N + l ro - 1 , -N + l m t| = a -2(A r -l) c -2( j V-l) g - 2(j V-l)(,+ m -l) 
, • • -, -N + l m -\ -N + 2 mt | = ^W-%-2(N-l) q -2(N-l)(v +m -l 



| -j^m — 1 2 m_ 

| -ym — 1 2 m — 

Difference Formula IV 



v.m 
'N ■ 



|0„ +1 , Of) ■ 



• ,N = 2„| = C-^ N - 1 \alq 21 ' - l)-(N-l) q -2(m-l)(N-l) ^+ 



+l,m 



0„, • ■ ■ ,N = 2„| = c-^-^agg 2 " - i)-(Ar-D g -2(m-i)(i V -i) ^+1,™ 

Difference Formula V 

|U„ +1) 0„, • ■ ■ ,N = 2„| = (1 - a 2 ^)-^" 1 ) V^ +1 ' m , 
IWi, 0„, • ■ ■ ,N = 2„| = (1 - alq^y^-Vafa 2 " V^' m ■ 



Table 1: Data for the proof of difference formulas 



Difference Formula 


Contiguity Relation 


Determinant 


Difference Formula I (97 


) 


G™ + Y = G™ +1 + a 2 q 2 ^+ m h 2 G^ 


eq.(|2 




Difference Formula 11(98) 


\J V — — q C l_z„ 


eq-(l 




Difference Formula 111(99) 


rim sim „2„2„2fj/+m— 2)/--rm— 1 
~ ^v-X — a c 1 W-l 


cq.(|3 




Difference Formula IV(lOOj) 


r-mi+1 rim _ r 3. r ,1m(\ _ rl 2v\/~tm 


eq.@ 


Difference Formula V(1CT 


) 


Lr„+l - aoq = (1 - a q )G„ 


eq.fl) 



Difference formula I has been proved in Lemma 2.5. Other formulas are proved in a similar manner. Table | 
shows the the contiguity relation and determinant to be used for the derivation of each difference formula. 



Next we give the list of Plucker relations which are necessary for the proof of Proposition 2.4 
Pliicker Relation I 







|*i- 1 

|*i, o 



m+1/ -i m 



x 10 



m i m 



m i m 



x 1 



m+1/ -i m 



I m+1/ r\m 



..,N-l m | x 1™,. 



•,N-l ro , N r ' 
,N- l m , N" 
,N-l m , N r ' 



Plucker Relation II 



= |*2, 1 
1*2, 



m+1/ -i m 



•••,N 



x |0 m , 1™ 



m i m 



..,N- 



x 1 



m+1/ -i m 



Pliicker Relation III 

= 

+ 

Plucker Relation IV 



I m+1/ r\m 



I m -i m+1 



•,N 



x |* 2 , 1"Y 



,N- l 1 
,N-1" 

,n- v 



N' 

N" 



|*2, 1"\ l r 
U, m+1 , 



m+1 i m+1 



o 



|*2, *1, l m+1 , 
|*2, I", l m+1 , 
1*1, 1™, l m+1 , 



,N-1 
.-,N- 
••,N 



,N- 

,N- 



m+1 I 



X 



m+1 -i m+1 



I m+1 I 



I m+1 1 



x \r 

X 1*2, 1 



I m+1 



N _!"•+! N ™+1| 



m+1 



,N 
,N 



I m+1 



I m+1 



N 



m+1 1 



I m+1 1 



I m+1 I 



I m+1 I 



X |l r 

X |(/3 1; 1 
X 1*2, 1 



I m+1 



m+1 



m+1 



,N 
,N 
,N 



I m+1 



I m+1 



I m+1 



m+1 1 



m+1 1 



Pliicker Relation V 



= \<p 2 , <p x , m ,- 
- |*2, m+v , 0' 

+ U, o m+v , o r 



,N-2 m | x |0 ro+1/ , 0" 
••■,N-2 m | x fa, 1 
•••,N-2 m x \<p 2 , r 



Plucker Relation VI 



= \o m ~\ 


~i 171— 1 
1 


•,-NH 


- 1 


- |0 m "\ 


-i m — 1 

5 


•,-NH 


- 1 


+ \o m ~\ 


-i m — 1 

5 


■,-NH 


-1 



m— 1 
m — 1 



,N-2 m , N-l m | 
,N-2 m , N-l m | 
,N-2"\ N-l m | 

N+1 m-l _ N 



, m+lf 



(102) 



(103) 



(104) 



(105) 



(106) 



*i| 



N + l mT x |-1 

m— 1 i „ i im— 1 



*l| x 



m — 1 

N- 



-N 



-N- 



■ mfi 



(107) 



Pliicker Relation VII 



Plucker Relation VIII 



- \ipi, U , l„,---,N-lJ x 



•,N-1„, NJ 



■,N-1„, N„ 



x lipx, l u , ■ ■ ■ ,N - 1„, NJ 



(108) 



= \(f2, (pi, lv , • • • ,N - l v \ x lOv+i, 1„, 



,N-1„, N„| 



|V? 2 , O^+i, 1 W ,---,N-1 I/ | X |</?i, l I/ ,---,N-l„,N I/ | 
O^+i, ,N - 1„| x |</?2, T„, • • • ,N - 1 V ,N V \ 



(109) 



Table 2: Data for the proof of bilinear equations 



Bilinear Equation 


Plucker Relation 


Difference Formula 


eq.( 


H) 


Pliicker Relation I (: 


02) 


Difference Formula I (97) 


eqs.© and @) 


Pliicker Relation II ( 
Pliicker Relation III 


103) 
(104) 


Difference Formula I (|97|) 
Difference Formula II (|98|) 


eq.( 


H) 


Pliicker Relation IV ( 


105) 


Difference Formula II (ps|) 


eq.( 


H) 


Pliicker Relation V( '. 


.06) 


Difference Formula I (197) 


eq.( 


H) 


Plucker Relation VlfllQ7|) 


Difference Formula III (p9) 


eq.( 


H) 


Plucker Relation VII(|l08|) 


Difference Formula IV (^00|) 


eq.( 


H) 


Plucker Relation VIII([l09|) 


Difference Formula V (|l()l|) 



Here, (pi (i = 1, 2) are any column vectors which we specialize as 









f ) 









<P1 = 




, (f2 = 















V o J 







(110) 



in the proof of bilinear equations. 

In Section 2.2.2, we have shown how eq.( f42|) ca n be derived from the Plucker relation I and difference formula 
I. The other bilinear equations in Proposition |2.4| are obtained similarly by using the data described in Table 0. 
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